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Abstract 

We show that in the framework of gravity and in the linearized 
approach it is possible to obtain spherically symmetric stationary states 
that can be used as a model for galaxies. Such approach could represent 
a solution to the Dark Matter Problem. In fact, in the model, the Ricci 
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curvature generates a high energy term that can in principle be identi- 
fied as the dark matter field making up the galaxy. The model can also 
help to have a better understanding on the theoretical basis of Einstein- 
Vlasov systems. Specifically, we discuss, in the linearized gravity, the 
solutions of a Klein- Gordon equation for the spacetime curvature. Such 
solutions describe high energy scalarons, a field that in the context of 
galactic dynamics can be interpreted like the no-light-emitting galactic 
component. That is, these particles can be figured out like wave-packets 
showing stationary solutions in the Einstein- Vlasov system. In such ap- 
proximation, the energy of the particles can be thought of as the galactic 
dark matter component that guarantees the galaxy equilibrium. Thus, 
because of the high energy of such particles the coupling constant of the 
i?^-term in the gravitational action comes to be very small with respect 
to the linear term R. In this way, the deviation from standard General 
Relativity is very weak, and in principle the theory could pass the Solar 
System tests. As pertinent to the issue under analysis in this paper, we 
present an analysis on the gravitational lensing phenomena within this 
framework. 

Although the main goal of this paper is to give a potential solution 
to the Dark Matter Problem within galaxies, we add a Section where we 
show that an important property of the Bullet Cluster can in principle be 
explained in the scenario introduced in this work. 

To the end, we discuss the generic prospective to give rise to the Dark 
Matter component of most galaxies within extended gravity. 

PACS numbers: 95.35. +d, 04.50.Kd. 

Keywords: galactic high energy scalarons; Einstein- Vlasov system; Dark 
Matter. 

1 Introduction 

The accelerated expansion of the universe, that is today observed, suggests that 
cosmological dynamics is dominated by the so-called Dark Energy field which 
provides a large negative pressure. This is the standard picture, in which such 
new ingredient is considered as a source of the right hand side of the field equa- 
tions. It should be some form of non-clustered non-zero vacuum energy which, 
together with the clustered Dark Matter, drives the global dynamics. This is 
the so-called ^^concordance model" (ACDM) which gives, in agreement with the 
Cosmic Microwave Background Radiation (CMBR), dim Lyman Limit Systems 
(LLS) and type la supernovae (SNela) data, a good framework to understand 
the today observed Universe. However, it presents several shortcomings as the 
well known "coincidence" and "cosmological constant" problems [IJ. An alter- 
native approach is to change the left hand side of the field equations, and check 
if observed cosmic dynamics can be achieved by extending general relativity 
131 m [51 [5] . In this different context, it is not required to search candidates for 
Dark Energy and Dark Matter, which till now have not been found. Rather, one 
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can only stand on the "observed" ingredients: curvature and baryon matter, to 
account for the observations. Considering this point of view, one can think of 
that gravity is not scale-invariant. Such an assumption opens a room for alter- 
native theories to be introduced [3 |H1 13 • In principle, the most popular Dark 
Energy and Dark Matter models can be achieved by considering f{R) theories 
of gravity where R is the Ricci curvature scalar. 

In this picture, even the sensitive detectors for gravitational waves like bars 
and interferometers (i.e. those which are currently in operation and the ones 
which are in a phase of planning and proposal stages) [T^ could, in principle, be 
important to confirm or rule out the physical consistency of general relativity or 
of any other theory of gravitation. This is because, in the context of Extended 
Theories of Gravity, some differences between General Relativity and the alter- 
native theories can be pointed out as far as the linearized theory of gravity is 
concerned pT[ IT^ . 

In the general picture of high order theories of gravity, recently the 
theory, which was originally proposed by Starobinski [T3], has been analysed 
in various interesting frameworks, see [TH [THl UHl E] for example. Specifically, 
the non-singular behaviour of this class of models is discussed in [M]. In [15] 
R^ inflation is combined with the Dark Energy stage and in [TB] an oscillating 
Universe, which is well tuned with some cosmological observations is discussed. 
Finally, in [T7] the production and potential detection of gravity-waves from 
this particular theory has been shown. 

It is also quite important to emphasize that the R^ is the simplest one 
among the class of viable models with i?™ terms in addition to the Einstein- 
Hilbert theory. In Ref. [F, it has been shown that such models may lead to 
the (cosmological constant or quintessence) acceleration of the universe as well 
as an early time era of inflation. Moreover, they seem to pass the Solar System 
tests, i.e. they have the acceptable newtonian limit, no instabilities and no 
Brans-Dicke problem (decoupling of the scalar) in the scalar-tensor version. 

In this paper it is shown that in the framework of R^ gravity, in the linearized 
approach, it is possible to obtain spherically symmetric and stationary galaxy 
states which can be interpreted like an approximated solution to the Dark Mat- 
ter Problem. In fact, in the proposed model the Ricci curvature scalar generates 
an energy term that can in principle be identified as the Dark Matter content 
in the galaxy. The model can also help to have a better understanding of the 
physics of Einstein- Vlasov systems pS] . 

As the observed gravitational lensing phenomena are well related to the S- 
type galaxies in the Hubble sequence, here we briefly discuss in Section 4 the 
implications of this theory for observations of bending of light by foreground 
galaxies. 

The main goal of this paper is to give a potential solution to the Dark Matter 
Problem within galaxies. In other words, we do not ban the presence of a Dark 
Matter component in our Universe which can explain properties of clusters of 
galaxies which cannot be achieved by alternative gravity theories. In any case, 
we add a Section where we show that an important property of the bullet cluster 
can be, in principle, explained in the scenario given in this paper |44] . 
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At the end of the paper, we discuss the general possibihty to give rise to a 
Dark Matter within extended gravity. 



2 The field equations and the hnearized approach 

Let us consider the high order action (TSl [111 ITT] 

S = J d^X^{R + bR^+Cm). (1) 

In order to avoid confusion, in ah this paper b represents the couphng constant 
of the term. 

The action ([1]) is a particular choice with respect to the well known canonical 
action characterizing General Relativity (the Einstein - Hilbert action [IS]), 
which reads 

f d\^{R + C^). (2) 



If the gravitational Lagrangian is not linear in the curvature invariants, the 
Einstein field equations has an order higher than second. For this reason, such 
theories are often called higher-order gravitational theories This is ex- 

actly the case of the action ([1]). Note that in this paper we work with SttG = 1, 
c = 1 and h — \^ while the sign conventions for the line element, which generate 
the sign conventions for the Riemann/Ricci tensors, are (— , -|-, -|-). 

By varying the action ([T|) with respect to g^^, (see refs. [TTJ [TCI [TT] for a 
parallel computation) the field equations arc obtained: 

G^, + b{2R[R^, - + 

-2R,^,, + 2g^,UR} = T^™' 
with the associated Klein - Gordon equation for the Ricci curvature scalar 

DR^ E^{R + T), (4) 

which is obtained by taking the trace of equation ([3]), where □ is the d'Alembertian 
operator and the energy term, E, has been introduced for dimensional motiva- 
tions: 

thus, b has to be positive |16) . 

In the above equations is the standard stress-energy tensor of the mat- 
ter. Note that General Relativity is obtained for 6 = in Eq. (3). 

Because we want to study interactions between stars at galactic scales, the 
linearized theory in vacuum (T^l^^ = 0), which gives a better approximation 
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than Newtonian theory, can be analysed, by considering a small perturbation of 
the background, which is assumed to be given by a Minkowskian background. 

In the linearization procedure we will perform a computation very similar to 
the ones in |16[|17) . but with some difference that will emphasize the importance 
of the Ricci scalar for the purpose of this paper: to provide a curvature-inspired 
model for galaxies. 

Putting 

g^Li' = Vf^iy + (6) 

to first order in h^i,, calling R^^pa, R^w and R, respectively, the linearized 
quantity which correspond to R^i^pa, R^i, and R, the linearized field equations 
are obtained as [111 HI] '■ 

Rfii^ - ^Tlpu = -df^d^bR + rif^yUhR 

UR = E'^R. 

Since R^^pa and Eqs. ([7]) are invariants for gauge transformations of the 

kind [ig Hz] 

(8) 

bR bR' = bR; 
then the gauge transformation 

hp^ = - ^rjp^ + Tjf,„bR (9) 

can be defined. Thus, by considering the transformation rule for the param- 
eter 

De,=d^hp,, (10) 

a gauge analogous to the Lorenz one for electromagnetic waves can be chosen, 
which then reads 

a'^ V = 0- (11) 

In this way, the field equations then read like 

nJip, = 0, (12) 

DbR = E^bR. (13) 
Solutions of eqs. and ([T^ are plan waves: 

hf,u = Apt.i'f) expiip/^xp) + c.c. (14) 
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bR ~ a{jf') e:iqy{iq'^xp) + c.c. (15) 



where 



P' 



{uj,'f), uj=p=\'f\ 



(16) 



In Eqs. (|12p and (|T4|) the equation and the solution for the tensor waves, 
exactly like in general relativity, have been obtained [121 HZ], while Eqs. (|13p 
and (1151) are, respectively, the equation and the solution for the third mode 
stemming from the curvature (see also [T6l [T7] V 

The fact that the dispersion law for the modes of the field bR is not linear 
has to be emphasized. The velocity of every ^^ordinary" (i.e. which arises from 
General Relativity) mode hfj,^, is the light speed c, but the dispersion law (the 
second of Eqs. (|16l) ) for the modes of bR is that of a wave-packet [16l[T7]. Also, 
the group-velocity of a wave-packet of bR centred in is 

= -. (17) 
From the second of Eqs. (ITBl) and Eq. ([T7]) it is simple to obtain: 



Then, it is also written like [161 HZI 



VG = . (18) 



E 



V/(1 - 4)^- (19) 



Now, the analysis can remain in the Lorenz gauge with transformations of 
the type Dei, = 0. This gauge gives a condition of transverse effect for the 
ordinary part of the field, i.e. k^^A^i, = 0, but does not give the transverse 
effect for the total field h^^. 

From Eq. ^ it reads 

- h ~ 
h^^y = h^^y - 2 ^M"^ + ViJ^ubR. (20) 

At this point, if one were still to remain within general relativity one would 
impose [25] . 

Ue^' = 

(21) 

which would give the total transverse effect of the field. However, for the 
present case this is impossible. In fact, applying the d'Alembertian operator to 
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the second of Eqs. (PT|) and using the field equations ([T^ and it comes 

out that 



9^ne^ = (22) 

which is in contrast with the first of Eqs. (pij) . In the same way, it is possible 
to show that no linear relation exists at all between the tensor field h^i, and 
the linearized term hR stemming from the curvature scalar. Thus, a gauge in 
which hfj^i, is purely spatial cannot be chosen (i.e. one cannot put /i^o = 0, see 
Eq. (EOl))- 

Nonetheless, the traceless condition to the field h^^ can be written in the 
form: 

□e'' = 

(23) 

These equations imply 

5'^ V = 0- (24) 
To save the conditions d^Ji^'^ — and h — {)^ transformations like 



a^De^ = 



^^€^' = 



(25) 



can be used and, by taking "p" in the z direction, a gauge in which only An, 
^22, and Ai2 = A21 are different from zero can be chosen. The condition h = 
gives All = -^22- 

Now, by substituting these equations in Eq. (EU)) . one gets 

h^^^t, z) - A+{t - z)e'^+J + (t - z)e%^ + hR{t, z)r,,,,. (26) 

The term A+(t — z)e^,t^ + A'^{t — z)e{i^^ describes the two standard polar- 
izations of gravitational waves obtained from general relativity, while the term 
hR{t, z)?7^y represents a mode arising from the curvature term associate to the 
i?2 high order theory pBfTT] . 

In other words, in the R? theory of gravity, the linearized Ricci scalar is a 
source of a third polarization mode for gravitational waves which is not present 
in standard general relativity. This third mode is associated to a "curvature^^ 
energy E (see Eq. (jU). 

We also recall that the idea of considering the Ricci scalar as an effective 
scalar field [scalaron] arose from Starobinski |13] . 
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3 Application to the Einstein-Vlasov system: sta- 
tionary galaxies 

Now, we are going to discuss a model of stationary, spherically symmetric galaxy, 
assuming that the dynamics of the matter, i.e. of the stars making of the galaxy, 
is described by the Enstein-Vlasov system [TB]. In this way, the gravitational 
forces between the particles of the system, i.e., a galaxy, will be mediated by 
the third mode of Eq. (l26l) . i.e. by the spacetime curvature. Thus, the key 
assumption is that in a cosmological context such a mode, which is given by 
the (linearized) spacetime curvature, becomes dominant at galactic and cosmo- 
logical scales (i.e. A^,A~ <C bR). In this way the "curvature" energy E can 
be identified as the dark matter content of a galaxy of typical mass-energy: 
E ~ 10*^(7 ^5\, in ordinary c.g.s. units. These two assumptions constitute 
an adaptation to modelling a galaxy of the assumptions related to the model 
describing an oscillating Universe in |16) . 

Now, let us emphasize an important point. Assuming E ~ lO^^g, from 
Eq. © we get b ~ 10"^'* cm^ in natural units. Thus, in our assumption, the 
constant coupling of the term in the gravitational action results infinitesimal 
with respect to the linear term R. In this way, the variation from standard 
General Relativity is very weak, thus the theory can pass the Solar System 
tests. Regarding this important issue, it is important to provide citations of 
precedent work illustrating this and explicitly show that the bounds there agree 
with the preferred values of the number E. The key point is that as the effective 
scalar field arising from curvature is very energetic, then the constant coupling 
of the the R^ nonlinear term — )■ [31]. In this case, the Ricci curvature, which 
is an extra dynamical quantity in the metric formalism, must have a range 
longer than the size of the Solar System. An important work is ref. |32| . where 
it is shown that this is correct if the effective length of the scalar field I is 
much shorter than the value of 0.2 mm. In such a case, the presence of this 
effective scalar is hidden from Solar System and terrestrial experiments. The 
value of E that we are assuming here guarantees the condition I ^ 0.2 mm. 
Another important test concerns the deflection of light by the Sun. This effect 
was studied in R^ gravity by calculating the Feynman amplitudes for photon 
scattering, and it was found that, to linearized order, this deflection is the same 
as in standard General Relativity 3^. 

On the other hand, Eq. (|15p guarantees that spacetime curvature does 
not remain too small and thus it can, in principle, become the Dark Matter 
component in the galaxy 24 . 

The model that we shall discuss is similar to the one introduced by Nord- 
strom in [21 . The relevance of the Einstein-Vlasov model was emphasized in 
the 1992 in a famous paper by Rein and Rendall [18]. And the following results 
will be obtained adapting the ideas introduced in [TS] |22l |24j |26l [28] . 

In the hypothesis ,A^ <C 6i?, the spacetime of our model will be given 
by the conformally flat metric 
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ds^ = [1 + bR{t, z)]{dx^ + dy^ + dz^ - dt^). 



(27) 



Note: in general, conformal transformations are performed by rescaling the 
line- element like [53] 

3a/3 = e*5r„^. (28) 
Here we choose the scalar field as being 

^ = bR (29) 

which also implies 

e* = 1 + bR, (30) 

in our linearized approach. Thus, the conformal transformation which trans- 
lates the analysis into the conformal frame (Einstein frame, see 0) is performed 
by the spacetime curvature. The potential of using such conformally flat line- 
element in the perspective of describing an oscillating Universe has been dis- 
cussed in [TB] . 

The condition that the particles in the spacetime make up an ensemble with 
no collisions (P7)) is satisfied if the particle density is a solution of the Vlasov 
equation [HI [H Hi HH] 

dtf+j,d.^f~Tl,i^d^.f = Q, (31) 

where are the ChristofFel coefficients, / is the particle density and is 
given by p°(a = 1, 2, 3) according to the relation [181 122< i26, ,28) 

9^..P''P'' = -1- (32) 

Eq. (I32p implies that the four momentum lies on the mass-shell of the 
spacetime (greek indices run from to 3) [TSJ I^Hl HH] ■ 

We recall that, in general, the Vlasov-Poisson system is given by [l8l[22l[26l 

[28] 

dtf + v- sj^f - \J^U ■ \7„f = 

\j .U ^4np (33) 

p{t,x) = J dvf{t,x,v), 

where t denotes the time and x and v the position and the velocity of the 
stars. The function U = U{t, x) is the average Newtonian potential generated 
by the stars. This system represents the non-relativistic kinetic model for an 
ensemble of particles with no collisions interacting through gravitational forces 
which they generate collectively [l8l [22l [26 l [28] . 
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Thus, such a system can be used for a description of the motion of the stars 
within a galaxy, if stars are considered as pointhke particles, and the relativistic 
effects are negligible [TFl [551 I^Hl [13 • In this approach, the function f{t,x,v) in 
the Vlasov-Poisson system is non-negative and gives the density on phase 
space of the stars within the galaxy. 

In other words, we are going to discuss the solutions of a Klein- Gordon 
equation for the spacetime curvature like galactic high energy scalarons, i.e. 
particles that can be figured out like wave-packets, and stationary solutions in 
terms of an Einstein- Vlasov system will be shown. 

In such approximation, the energy of the particle will be seen like the Dark 
Matter component that guarantees the galaxy's equilibrium. It is stressed that 
such approximation is not as precise as we would, but it could be a starting 
point for further robust analysis. 

The Vlasov equation (|3ip implies that the function / is constant on the 
geodesic flow over the spacetime (1571) . The ChristofFel coefficients of such a 
spacetime are obtained from (note that as we are working in the linearized 
approach, in the following computations only terms up to first order in the 
linearized curvature R will be considered while high-order terms will be assumed 
equal to zero) 

r^. = Us^d^bR + S'^^d.bR - -J—^g^.d'^bR). (34) 

In this way, the Vlasov equation in the spacetime defined by the line-element 
(pTl) becomes 

p" 1 
dtf + -r^d^af - - 
2 

Now, let us recall that two quantities are important for the Vlasov equation 
in a curved spacetime [T51 [5^1 HH] • The first is the current density 

N^ = -j^V9P'f (36) 
and the second is the stress-energy tensor 

T^- ^- J ^y/^p'^p-'f. (37) 

Here g is the usual determinant of the metric tensor, which in the case of 
the line-element (|27p is given by 

.9 = 1 + AbR. (38) 

Keep in mind that both N'^ and T^"^ are divergence free (conservation of 
energy): 



2ip^d^bRy- 



d^bR 
(1 + 2bR)p^ 



(35) 
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V^TV^ = 0, 

(39) 



S/f.Tf" = 0. 

The mass shell conditions (|32]) can be rewritten as 



/ = V (1 + bR)-^ + SabP'^p''. (40) 

From the ChristofFel connections computing the Riemann tensor, Ricci 
tensor and Ricci scalar, the "effective"^ Einstein field equations 

G^. = T^., (41) 
can be obtained together with the "effective''' Klein - Gordon equation 

UhR ^ -T, (42) 
where T = is the trace of the stress-energy tensor. 

To simplify the computations the analysis can be performed in a conformal 
frame. Thus, rescaling the stress-energy tensor in the form 

= (1 + 36^)^^^ (43) 

one obtains 

(1 + 'SbR)T. (44) 

Then, equation (H^ becomes 

DbR = -n. (45) 

It should be noticed that the particle density is still defined on the mass shell 
of the starting line-element ga^. In order to remove even this last connection 
with the starting frame (Jordan Frame) we can rescale the momentum as 

p'.={l+^-§)p' (46) 
and we can define the particle density in the Einstein frame as 

Mt,x,p,)^ f(^t,x,{l-^)p}j . (47) 

Hence, we can write our adaptation of the Vlasov system in the following 
form 

DbR={l + 2bR) f ^Mt,x,p,), (48) 
J Pi 



p" = \ I + SabPtpl- (49) 
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dtf* + 4^-"/* - ^[p'id^.bRpt + d%R]dp.U - 0. (50) 

Because we want to restrict ourselves to spherical symmetry in the present 
approximation, the line-element (|27p can be rewritten as 

ds^ = [l + bR{t,r)]{dr^ -df). (51) 

In this equation r is the radial coordinate. Thus, in spherical coordinates, 
equations (|48l) . (|49|) and ([50l) can be written as 



d'^bR I d f d ~ y\ „ ~, , , 

/i(t,r)= / -^=^fit,x,p), (53) 




= 0, 



(54) 

where the suffix * has been removed for the sake of simplicity, and we have 
denoted by p the vector p — (pi,p2,P3) with p^ = and also defined x for 
the vector Xi = {xi,X2,X3). 

Then, in searching for stationary states, and following [TB], we can call A 
the wavelength of the "galactic" gravitational wave (|26p . i.e. the characteristic 
length of our gravitational perturbation, and assume that A ^ d, where d is the 
galactic scale, i.e. d ~ 10^ light-years [25]. In this way, the gravitational wave 
is "frozen-in" with respect to the galactic scale. 

Thus, the system of equations which defines the stationary solutions of eqs. 
([5^ . ([55)) and (IMl) . for our galaxy model is 

^ ^ '^-bRrA ^{l + 2bR)ii{r), (55) 



dr \dr 

^,{r) = / ^^/(x,p), (56) 

p-dj --^bR[{p-x)p + x]-dj,f = 0, (57) 
r dr 

In other words, the idea is that the spin-zero degree of freedom introduced via 
the addition of an R^ term in the gravitational Lagrangian may be a candidate 
for the dark matter. This new degree of freedom is termed the scalaron |13) . 

Thus, we assume that within a galaxy the dominant contribution to the 
curvature comes from the new degree of freedom (as it is also the case to first 
approximation in the general relativity -f cold dark matter case). We deduce 
this contribution via the scalaron field equation which itself has a baryonic 
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source term. The baryons themselves are then taken to evolve according to a 
coUisionless Bohzmann equation, propagating on a background perturbed by 
the scalaron. The coUective set of equations can in principle be solved for given 
initial data. 



4 The gravitational lensing 

The attentive reader notices that, in principle, as the metric is conformally 
flat, this could mean that the scalaron Dark Matter in galaxies will, in and of 
itself, result in no gravitational lensing of light rays. This could be a serious 
problem for the model. This is because there exists substantial evidence for 
gravitational lensing by the large scale structure of a galaxy beyond that due to 
the constituent baryons alone. For instance, the observed absence of lensing by 
the effective dark matter in a class of metric formulations of Modified Newtonian 
Dynamics formed the basis for a tentative no-go theorem for such theories [30] . 

Actually, we show that in the proposed model the gravitational lensing can 
be, in principle, obtained like an effect of spacetime curvature, i.e. due by the 
Ricci curvature scalar. 

In our linearized approach, gravitational lensing can be described in a lo- 
cal Lorentz frame perturbed by the first order post-Newtonian potential j27j . 
Calling V such a potential one can define a refractive index [57] 

n=l-|-2|y|. (58) 

Some clarifications are needed concerning this point. In the usual Geomet- 
rical Optics, the condition n > 1 implies that the light in a medium is slower 
than in vacuum [33]. Then, the effective speed of light in a gravitational field 
is expressed by [36] 

2\V\. (59) 

n 

Thus, one can obtain the Shapiro delay [37^ by integrating over the optical 
path between the source and the observer: 

observer 

2\V\dl. (60) 

ource 

The situation is analogous to the prism [36] . 

Now, we recall that, in the weak field approximation, the connection be- 
tween the post-Newtonian potential and the linearized theory is given by the 
goQ component of the line-element [TH] : 

goo « 1 + 2V. (61) 
Thus, from eq. (I?71) one obtains 

(62) 
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and 



n^l + 2\bR{t,z)~l\ (63) 

is the equation which can, in principle, be used to discuss the gravitational 
lensing in our model. Then, in this case, the gravitational lensing is performed 
directly by spacetime curvature, i.e. by the Ricci scalar. 

For a sake of completeness, we discuss the gravitational lensing from another 
point of view. 

The condition A ^ d, on the characteristic length of the gravitational per- 
turbation, guarantees that the Ricci curvature remains frozen, i.e. constant,^^ 
with respect to the galactic characteristic distance scale. Thus, we put 

bR = K{a constant!). (64) 

We can search deviations from standard General Relativity within the galaxy 
by considering a spherically symmetric Schwarzschild-like metric generated by 
the ordinary (barion) galaxy mass M with the corrections that are generated 
by curvature [38] : 

ds'^ = - exp{-\r)df + exp(Ar)dr2 + r^{de^ + sin^ Odip^) (65) 
Then, following [35], the Ricci scalar is given by: 

R = exp{~\r){X" + + (66) 

where ' stands for the derivative with respect to r. 
Putting the condition ([64)) in Eq. ([66[) one gets [38] : 

and, by choosing a = —2M, (3 — Oin analogy with the standard Schwarzschild 
metric, one gets [111 13H] 



ds^ = -(i-f^-^r^)dt^ + {l^ — + -^r^)-'dr^+r^{de^+sin^9d^^). (68) 
r 12b r 126 

The physical interpretation implies that within some radius Tmin deviation 
from General Relativity are highly suppressed and we say that the object is 
screened [39[ 1401 141] . This is because we are assuming that, within the radius 
'"min, the spacetime curvature which is due by the ordinary, i.e. barion, mass of 
the galaxy dominates with respect to the intrinsic constant curvature R = K . 
On the other hand, for values r > r-n-nn, the spacetime curvature which is due by 
the intrinsic curvature dominates with respect to the spacetime curvature due 
by the barion mass of the galaxy. 

In order to clarify the issue of the bending of light in this R^ scenario, in what 
follows, and as matter of illustration, we will reproduce a piece of the discussion 
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presented by Smith in Ref. |41| . which, in comparison with our discussion at 
the beginning of Section II, appears as a generalization of that digression (see 
the Appendix). 

Thus, for a general action S = J d^x^^—g{R + f{R) + Cm) screening occurs 
within a radius rmin implicitly given by |39l |40l |41] 



|/'(i?max)| <p(r„,in)rLn, (69) 

where p is the local value of the density. 

In our case it is f'{R) = 2bR and i?max = -y- Then, Eq. ([69]) becomes 

K < ^/3(r„,i„)rii„. (70) 

Now, following [41] , we can use Eq. (|70l) to further discuss the gravitational 
lensing. 

Let us consider a galactic Navarro, Frank and White (NEW) halo of the 
form gni E] 

p{r)^PcSc (71) 
r(r + rg) 

where is the scale radius and pc the critical density of the universe [301 E] ■ 
By assuming the mass-concentration relation [JI] |33J 23] 

^ (72) 

.12 X 10i2/i-iM®y ' ^ ' 

where z is the redshift, M@ the mass of the Sun and h is the Hubble param- 
eter in units of lOOfcm/ (sMpc), the amplitude Sc, that relates the concentration 
to the virial radius with an overdensity A = 119 is given by [TO H5] 

A 

'^^^3[ln(l + c)-(^)]- 

Eq. ([72]) reduces the NEW halo profile to a one-parameter family which is 
taken to be dependent on the virial mass, M [IT]. Eor r < rg pnfw scales like 
r^^ , while for r > Pnfw scales like r~'* [41| . Thus, the innermost point at 
which deviations from general relativity are suppressed in will occur at the scale 
radius rg ■ Considering Eq. (|70p togheter with the NEW density profile on sees 
that halos with masses which satisfy 



K > pMM) (74) 

As 5c{M) decreases with decreasing M [41], an upper limit is settcd to 
the mass of halos which can be screened given i.e. given the value of the 
Ricci curvature, see Eq. ([64| . Eor masses below this threshold the theory 
strongly deviates from general relativity and strong lensing around these halos 
is permitted. We plan to further develop in this direction of research. 
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5 Can gravity explain the Bullet Cluster dy- 
namics? 



The scalaron has two aspects as a field (hke classical electric or magnetic field) 
and a particle (like quantized photon) . We treated the scalaron as a field in the 
Vlasov equation (see [H] for a recent cosmological application) whose scenario 
is in some sense similar to that of the MOND [JS]. In such a scenario, it seems 
to be difficult to explain the Bullet Cluster [46] although the Bullet Cluster 
could be naturally explained if the scalaron could be a heavy particle . 

Aside from this, we emphasize that the main goal of this paper is to give a 
potential solution to the Dark Matter Problem within galaxies. Therefore, we 
do not ban the presence of a Dark Matter component in our Universe which can 
explain properties of clusters of galaxies which cannot be achieved by current 
alternative gravity theories [46]. In any case, in this Section we show that an 
important property of the Bullet Cluster can in principle be explained in the 
scenario proposed in this paper. 

5.1 A note on the Bullet Cluster 1E0657-56 

Since 2006, the argument "But the Bullet Cluster ..." has come into scene in 
most discussions about the Dark Matter Problem in astrophysics. It relates to 
the observations of the post-collision clusters of galaxies in the source 1E0657-56 
(known as the Bullet Cluster) [IIT . It is seen that the X-ray emission is a mid the 
couple of colliding clusters, and thus is shifted with respect to the bulk of matter 
making up each group of galaxies, which resides in the outskirts. In most studies 
of X-ray emission from galaxy clusters the gas usually sits at its center. This 
offset can be understood by recalling that as the gas pass each other it interacts 
electromagnetically, heats up, ionize to emit X-rays, and hence is slowed-down 
during the collision. This may explain why the gas appears a bit behind in 
the heading direction of each of the group of galaxies involved. Meanwhile, 
non interacting Dark Matter components do only interact through gravity and 
thus pass each other unhindered without being slowed-down like the gas. In 
the image, the field purported to represent the vast of mass in each cluster, 
as compared to the X-ray emitting counterpart, is inferred from gravitational 
lensing analysis of images of background galaxies, under the assumption that 
general relativity is the theory of gravity (but be aware that most modified 
gravity theories describe correctly the phenomenon of gravitational lensing. See 
the Appendix). It follows the distribution of galaxies, not the gas. As the most 
massive component of the galaxy cluster is not centered at the gas, it is argued, 
upon that lens analysis, that the large part of matter in the cluster is close to, 
and around each galaxy group. And as the visible mass in the galaxies is not 
enough to account for the velocity dispersion of a galaxy cluster, it is suggested 
that there is Dark Matter, too. As noninteracting component, the distribution 
of the Dark Matter can pass through each other just like the galaxies, and that is 
why the bulk of mass should be found close to the galaxies in a cluster collision 
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like this. 

The analysis of the overimposed optical and X-ray images of the cluster can 
be interpreted as an endorsement on the existence of Dark Matter. This issue 
is based on a number of assumptions which are not totally accepted by the 
astrophysics community because the large set of problems for the Dark Matter 
hypothesis are, in any case, independent from the observations of the Bullet 
Cluster. In a different context, the collision velocity of the Bullet Cluster seems 
to be in disagreement with the standard concordance model of cosmology [IS] 
(the relative velocity of the clusters in the collision is an issue being addressed 
next in the framework of this gravity theory). At the same time, alternative 
theories of gravitation, while often said to fail in explaining dynamics of galaxy 
clusters, can account for them rather naturally. 

On the other hand, even if the Bullet Cluster can only be explained by 
invoking Dark Matter, the problems on small scales persist. The Bullet Cluster 
does not improve our understanding of the Local Group of Galaxies, nor on the 
dynamics of individual galaxies, (as is on focus in the present paper). The two 
arguments are absolutely independent. A conclusion could be that the lensing 
mass and the hot gas have a spatial shift, and hence the visible, hot gas cannot 
make up the bulk of mass in the system. In any case, whether the cause of this 
phenomenon is a nonstandard law of gravitation or missing mass is not that 
easy to unveil. 

5.2 Scalaron and the Bullet Cluster dynamics 

As pointed out above, it is a common opinion that alternative gravity theories 
cannot explain all the properties of the bullet cluster [A&^ , even if scalar- vector- 
tensor theory could, in principle, take into account the apparent discrepancy be- 
tween the gravitational lensing and X-ray maps of the colliding clusters, without 
requiring collisionless Dark Matter (471 [50] . 

Meanwhile, hydrodynamic simulations of the colliding galaxy clusters 1E0657- 
06 show that at a distance of 4.6 Mpc between the two clusters, an infall velocity 
of 3000 km/s is required in order to explain the observed X-ray brightness and 
morphology of the cluster [171 EH]- It has also been argued that such a high 
infall velocity is incompatible with predictions of the standard cosmology (the 
concordance ACDM model) [13 HI]. 

Next, we extend our scalaron model in order to show that the infall velocity 
of the two clusters can in principle be explained in the framework introduced 
in the present paper. To rigorously show that the scalaron gravity can cope 
with the Bullet Cluster dynamics we need first to modify our assumptions. In 
particular, we will assume that the characteristic length of the scalaron gravita- 
tional perturbation is much longer not only than the galactic scale, but than the 
distance between the colliding clusters too. In other words, the gravitational 
wave is "frozen-in" with respect to the Bullet Cluster's scale. In this way, we 
can consider the colliding clusters (the main cluster and the "bullet") like test 
particles in the gravitational field of the scalaron. 

In the linearized approach of this paper, the coordinate system in which the 
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space-time is locally flat has to be used and the distance between any two points 
(the colliding clusters) is given simply by the difference in their coordinates in 
the sense of Newtonian physics [12] ■ This frame is the proper reference frame 
of a local observer, which we assume to be located within the main cluster. 
In this frame gravitational signals manifest them-self by exerting tidal forces 
on the test masses. In other words, we assume that the space-time within the 
main cluster is locally flat with respect to the global curvature generated by 
the scalaron which is described by the line element (P7|) . By using the proper 
reference frame of a local observer the time coordinate xq is the proper time 
of the observer O and the spatial axes are centred in O. In the special case of 
zero acceleration and zero rotation the spatial coordinates the proper 

distances along the axes and the frame of the local observer reduces to a local 
Lorentz frame [19]. As the origin of the coordinates is located within the main 
cluster, the proper distances between the two clusters are the coordinates of the 
"bullet" . The line element is [H] 

ds^ = +{dx°f - Sdx'dx^ - 0{\x^^)dx'^dx^. (75) 

The effect of the gravitational force on test masses is described by the equa- 
tion 

x^ = -R}^k^x\ ^ (76) 

which is the equation for geodesic deviation in this frame [19, . i?ofeo 
linearized Riemann tensor and x'' are the coordinates of the "bullet" which 
represent the separation vector between the two test masses [IH] . 

To study the effect of the scalaron on the two clusters, i?Qj,o has to be com- 
puted in the proper reference frame of the main cluster. But, as the linearized 
Riemann tensor Rf^upa is invariant under gauge transformations [19', it can be 
directly computed from Eq. (|27p . 

From |19] we get 

Rfiuap — -^{dftd/^hau + dydah^j^p — dadph^„ — dfj^d^ha/s}, (77) 
that, in the case Eq. (I?71) . gives 

^070 = ^{d'^doRm-, + dod^RS^ - d'^d^Rrjoo - dod^^RS^}. (78) 

The different elements are (only the non zero ones will be written down 
explicitly) 

( dfR /or a = 7 = ] 

d'^doRT^o^ ^{ ^ \ (79) 

[ ~d,dtR for a = 3;7 = J 

( dfR for Q = 7 = \ 
dod^RS^ ^} \ (80) 

[ dtd.R for a = 0;7 = 3 J 
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—dfR for a = 7 = 

for a = 7 = 3 

—dfdzR for a 0; 7 = 3 

dzdtR for a = 3;7 = 



~dodaRd^= ~dfR for a = 7 
By inserting these results in Eq. (|78l) we obtain 



(81) 



(82) 



^010 



'.R 



??2 



(83) 



-^030 - ^{dlR- dfR). 

At the Bullet Cluster's scale we have homogeneity and isotropy, which imply 
dlR = 0. Therefore, Eqs. ^ become 



^010 



p2 



-^R 



(84) 



By using Eq. ([75)) . we obtain 



and 



(85) 
(86) 



z — —Rz. 
2 



(87) 



These are three perfectly symmetric oscillations of the scalaron (wave-packet). 
If one re- introduces the radial distance r, Eqs. (f85|) . (f86| and ([87| are summa- 
rized by 



-Rr. 



(88) 



Eq. 



can be solved with the perturbation method [19] 
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rit)^rmi + ^R{t)]. (89) 

Assuming that the wave-packet is in the contraction phase we get the relative 
infall velocity as 

v,^fau^j^r{t)^r{0)^Rit). (90) 
Therefore, in order to have consistence with the observations we need 

r{0)^R{t) - 3000 km/s. (91) 
Being r(0) 4.6 Mpc, the model will be consistent with observations if 

b ~, , 3000 km , „ „ i? i /„„x 

The assumption that the characteristic length of the scalaron gravitational 
perturbation is much longer than the distance between the colliding clusters 
implies that the frequency of the wave-packet has to be oje ^ lO^^^s^^. For 
example, for a value lue ~ 10~^^s~^ Eq. ((92|) guarantees theoretical consistence 
if 10"^ 

6 Final discussion and closing remarks 

It has been shown that in the framework of i?^gravity and in the linearized 
approach, it is possible to obtain spherically symmetric and stationary galaxy 
states which can be interpreted like an approximated solution of the Dark Matter 
problem. In fact, in the model the Ricci curvature generates a high energy 
term that can in principle be identified as the dark matter content of most the 
galaxies. The model can also help to gain a better understanding of the basics 
of the Einstein- Vlasov system. 

In other words, we have discussed, in the linearized i?^gravity, the solutions 
of a Klein-Gordon equation for the spacetime curvature like galactic high energy 
scalarons, i.e. particles that has been analysed like wave-packets, and we have 
shown stationary solutions in terms of the Einstein- Vlasov system. 

In this approximation, the energy of the galaxy-particle has been seen like 
the dark matter component which guarantees the equilibrium of galaxies. 

An important point is that, because of the high energy of such particle- 
galaxies, the constant coupling of the B? in the gravitational action comes out 
to be infinitesimal with respect to the linear term R. In this way, the deviation 
from standard general relativity is very weak, and, in principle, the theory could 
pass the Solar System tests. An analysis on the gravitational lensing has been 
also performed in Section 4. 

We recall that the main goal of this paper is to give a potential solution to 
the Dark Matter Problem within galaxies. Thus, we do not ban the presence 
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of a Dark Matter component in our Universe which can explain properties of 
clusters of galaxies which cannot be achieved by alternative gravity theories. 
In any case, in Section 5 we showed that an important property of the buUet 
cluster, i.e. an infall velocity of 3000 km/s between the two clusters, can be, in 
principle, explained in the scenario given in this paper. 
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Appendix: dynamics in f{R) theories of gravity 
and gravitational lensing 

We emphasize that hereafter we closely follow the paper by Smith [IT] , and the 
set of references quoted there. A generic branch of scalar-tensor theories in which 
the Brans-Dicke parameter uj = 0\ bring in unique predictions: they feature a 
peculiar effect known as the chameleon mechanism, by which the modifications 
to general relativity get quickly suppressed nearby sufficiently massive objects, 
such as galaxies. This rapid decay in the theory's gravitational effects present 
a unique size-dependent lensing footprint which is discussed next. 
The action for /(i?)-theories (functional of the Ricci scalar R) 



where Sm action for the matter fields. The field equations turn out to be 



In this class of theories the Ricci scalar becomes a dynamical quantity whose 
equation of motion is determined by the trace of the field equation. 



In the limit where / — >■ 0, Eq. ([M]) reduces to the standard relation R = —kT 
of general relativity. Here T is the trace of the stress-energy tensor. 

Using the trace equation to rewrite the gravitational equation of motion one 
obtains 




(93) 



[1 + /'(i?)] G^, + -g^ARf'iR)-f + 2Df'iR)] 
-V^V^/'(i?) = kT^,. 



(94) 




(95) 
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Teff = ^[^T + R] + ^ W{R) - /] • (96) 

Solutions to the Eq. ([Ml) determine the lensing predictions for the theory, 
and can be understood by rewriting the trace of the field equation in the form 



where 



dV 1 



dm-3V^^ + ^['-^'^^)l+'V- ^''^ 
Notice that for functions f{R) which reproduce the observed expansion his- 
tory of the universe, the minimum of this potential yields the general relativistic 
relationship between R and T, i.e, R = —nT. 

Hence, within a given object, a typical galaxy; for instance, far away from 
the center the scalar curvature starts off nearly at its asymptotic value, i?max5 
and evolves with radial distance from the center. If the object is too small 
compared to the wavelength of the gravitational perturbation, then R ^ i?max 
throughout the object, and deviations from general relativity will be of order 
unity. Meanwhile, if the object is large enough then the scalar curvature is 
forced to the minimum of its potential and R = —kT within some radius r^in- 
Within that radius deviations from general relativity are highly suppressed, and 
thus one says that the object is "screened''' . Finally, on the generic prospective 
for it to give rise to dark matter within extended gravity, we recall that a further 
argument in favor of the potential of modifying the gravity theory at the galactic 
scale so as to explain dark matter has been provided by Ellis [24] . 
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